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Spherical field theory is a new non-perturbative method for studying quan- 
tum field theories. It uses the spherical partial wave expansion to reduce 
CN ■ a general ci-dimensional Euclidean field theory into a set of coupled one- 

dimensional systems. The coupled one-dimensional systems are then con- 
verted to partial differential equations and solved numerically. We demon- 
strate the methods of spherical field theory by analyzing Euclidean <^> 4 theory 
in two dimensions. [PACS numbers: ll.10.Kk, 11.15Tk, 11.30Qc] 



oo ' 1 Overview 

We present a new non-perturbative method for studying quantum field the- 
ories. We describe a general procedure for converting problems in quantum 
field theory into initial- value problems described by partial differential equa- 
tions. Our derivation of the method starts with the functional integral for a d- 
dimensional Euclidean field theory. Making use of spherical 0((i)-symmetry, 
we decompose the fields as a linear combination of spherical partial waves. 
Re-interpreting the radial part of each partial wave as a distinct field in 
a new one-dimensional field theory, we rewrite the functional integral as a 
time-dependent Schrodinger equation, where the radial distance in the origi- 
nal theory serves as the time parameter. Although the resulting Schrodinger 
equation contains an infinite set of interacting partial waves, we find that 
high spin partial waves decouple if the original field theory is renormalized 
to remove ultraviolet divergences. Neglecting higher spin partial waves and 
estimating the corresponding error, we solve the remaining partial differential 
equation using standard methods. 
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The organization of this paper is as follows. We begin with a short discus- 
sion of one-dimensional systems and the connection between spherical field 
theory and coupled one-dimensional systems. Next we illustrate the method 
by analyzing normal-ordered Euclidean </> 4 theory in two dimensions. We 
derive the spherical Feynman rules for this theory, compute the two-loop 
self-energy, and compare the result with standard perturbation theory. We 
then consider non-perturbative values of the coupling constant and evalu- 
ate the self-energy numerically using finite-difference methods and diffusion 
Monte Carlo. Since our paper is intended as an introduction to the meth- 
ods of spherical field theory rather than a detailed analysis of Euclidean 4 
theory, we have not engaged in extensive numerical computations for our 
analysis here. Nevertheless, our rather modest calculations show clear evi- 
dence of spontaneously broken — > — </> reflection symmetry at the critical 
coupling value 

4 (1) 

where /j is the bare mass and A is defined by the quartic coupling ^(p 4 . Our 
value of the critical coupling is in close agreement with the results of a recent 
lattice study [1]. The methods we present in this paper are easily generalized 
to any (/-dimensional Euclidean field theory involving any number of scalar 
fields. The extension to fermionic and gauge systems is the subject of current 
work. 



2 Spherical fields 

We start with the Feynman-Kac formula, which connects the Schrodinger 
and path integral representations of (Euclideanized) quantum mechanics, 

(x F \Texp[-j^dtHj] \ XI ) (2) 



oc 



4>(ti)= x i 

4>(tF)=XF 



where 



Hj{t) = -\^ + V{q)-Jq (3) 
= H -Jq 
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and \xj) , \xp) are position eigenstates centered at xj, xp. The generating 
functional Z[J] is proportional to the ground state persistence amplitude, 



Z[.7]oc(0|Texp [- Too dtHj] |0> , 
where |0) is the ground state of H . For ti > t 2 , 



(4) 
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(5) 

The Feynman-Kac formula can be generalized to include more general time- 
dependent interactions, 



{x F \Texp dtH(t) 



oc 



cxp 
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(6) 



Ji)(!) 2 + mt) 



<t>(tl)=Xl 

4>{t F )=x F 



where 



(7) 



^(*) = -T^ + n«,*)- 

The central idea of spherical field theory is to treat a general ci-dimensional 
field theory as a set of coupled one-dimensional systems. To demonstrate, 
we consider Euclidean 4 theory in two dimensions, 



Z[J] = J £>0exp j-/d 2 x 



!£W) 2 + ^ 2 + !^ 4 -^ 

i 

Decomposing (j) and J into partial waves, we have 



(J)(t cos 9 , t sin 9) = 



E Ut)e ine , 



(8) 

(9) 
(10) 



J(tcose,tsm6) = J±- E Jn(t)e ine . 

V n=0,±l,... 

The decomposition of <fi into partial waves is a linear invertible transforma- 
tion. The Jacobian of the transformation (tempered with suitable ultraviolet 
and infrared regulation) is, 



det 



8<j>(t' cosfl/ sing) 

S4>„(t) 



det 



M( t - f ) e 
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This Jacobian is independent of (f> n (t), and we change integration variables 
to write Z\J\ as 



bl II Vfaexpt-ffdt 

n=0,±l,... 



n=0,±l,... 

"4!2 



t d<p- n d<j>n | /i 2 t 2 +n 2 , A, _ f T A 
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711 ,712 ,"3 

=0,±1,... 



(12) 

where 6 is an overall constant whose value is irrelevant for the purposes of 
our discussion. From (6) and (|7|), this functional integral corresponds with 
the Schrodinger time evolution generator 
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(13) 
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In order to define Z[j7"] in the Schrodinger language, we first restrict the 
t-integral in ( |12"1) to a finite interval, tj <t <tp. We then have 



z lti,t F ][J] oc n dx n da/„ ({<} n=0 ,±i, 

«/ n=0,±l,... 



where 



=o,±i, 



(14) 

{ X n} n =0,±l,~) • ( 15 ) 



We recover Z[J] by taking the limits tj — > 0, £p — > oo. The details of 
the quantum state at the boundary £p is rather unimportant since only the 
ground state projection survives as tp — > oo, similar to what happens in (4). 
In contrast, the go dependence of the quantum state at £/ = is significant. 1 
The proper initial state is 



\ a ) — / II d,X n \Xn\ n= Q±i 
J n=0,±l,... 

a state whose wavefunction is constant in g-space. 



(16) 



1 For n / the q n dependence is again unimportant and depends only on the ground 
state projection. 
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3 Correlators 

The generating functional, written as a power series, is 

where f{k 2 ) is the propagator, 

f(P) = fd 2 xe rk -*(0\<t)(x)m\0)- 



(17) 



(18) 



As a result of angular momentum conservation, the term quadratic in J can 
be decomposed as a sum over terms proportional to J- n J n , 



\ J]] Jdkdt x dtyt x t y k^— n (t x )^ n (ty)J\ n \(kt x )J\ n \(kty)f(k 



(19) 



n=0,±l, 



where Jj is the i th order Bessel function of the first kind. The two-point 
correlator for the n-th partial wave is 



(010^(^^(^)10)=^ ,^^ 



Z[0] ti<5J n (ti)i 2 5J r -„(« 2 ) 



J=0 



JdkkJ^kt^J^k^fik 2 ). 

(20) 

We will use this result in the next section when discussing perturbative spher- 
ical field theory. For the special case n = and t 2 = we have 

(O|0 o (ti)0o(O)|O) = JdkkJ (kt 1 )f(k 2 ) = J^e-^f(f) (21) 
= 2tt (O|0(r)0(O)|O), 

when |f| = t\. The two-point correlator is proportional to the two-point 0o 
correlator with one 0o centered at the origin. We can compute the two-point 
0o correlator using the Schrodinger formalism, 
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(22) 



More general correlation functions can be evaluated by partial wave expan- 
sion. For example the four-point correlator is given by 



ni,n 2 ,n 3 
=0,±1,... 



(oi 0(^)0(^)0(^)0(0) io) = ^ ifli g;r 3fl3) (oi 0„ 1 (ti)0n 2 (t 2 )0 ft3 (t 3 )0o(o) io) 

'■3 

(23) 
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where 

fj = (U cos 6i,ti sin 6i), i = 1,2,3. 
In the Schrodinger formalism, we can rewrite 

(o| <MO<M*2)<M*3)<fo(o) 1°) 

as (for t 1 >t 2 > t 3 ), 



(24) 
(25) 
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4 Perturbative spherical field theory 

For free field theory we have f(k 2 ) = p l + ^ 2 ■ From (20), 

(0| <t>- n {tl)<t>n{t2) |0) = ^(tl-t 2 )^M(^l)/H(^2)+^(t2-tl)^H(M2)/|n|(^l), 

(27) 

where and Ki are the i th order modified Bessel functions of the first and 
second kind respectively. In particular, we find 



(O|0oW0o(O)|O) = K o (^). 



(28) 



We can use (27) to generate the spherical version of Feynman diagrams. 
As an example, the two-loop process shown in Figure 1 gives us 



) 2 



0(t - t^Koifd^oifJth) 
+0(*i - QKoifrtJIoifrt) 



e{h - ^2)^1,2,3(^1)^1,2,3(^2) 

+6(t 2 - ^1)^1,2,3(^2)^1,2,3(^1) 

(29) 



where 



^1,2,3(3?) = h(x)I 2 (x)I 3 (x), 
K 1>2 , 3 (x) = K x {x)K 2 {x)K^x). 



(30) 



This diagram contributes to the two-point 0o correlator at order A 2 . Let 
A 2 ,j max (t) be the sum of all such 0(A 2 ) diagrams for the O correlator, with 
the restriction that each internal line correspond with some n such that 



K (iit 2 ), 
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\ n \ < 4ax' In view of relation (21) we expect that as J max ^ oo we should 
recover the regular two-loop perturbative result, 



1 /"l — OL\ 

A 2 / / rfa 2 rfai<5(ai+a2+Q3-l) 



967r 2 (fc 2 +/i 2 ) 2 J y aia2a :j k 2 +(aia2+a2Ce 3 +a 3 ai)fj, 2 ' 



lim ± fdtdete ukcose A 2 ,j_(t) 

(31) 

In Figure 2 we have plotted 



sj^(k 2 ) = f dtd6te ukcose A 2J _(t), J max = 0,l,2,3, (32) 

and 

s co(k ) «iQ2a3fc 2 +(«ia 2 +a2«3+a3ai)M 2 ' (33) 

in units where /i = 1. It appears that sj max (A: 2 ) does in fact converge towards 
Soc(k 2 )- We discuss this and related issues in the next section. 



5 High spin decoupling 

Partial waves with high spin correspond with high tangential momentum 



modes. We can see this explicitly in the coefficient of (p-n^n in (12]), 



(34) 

where the n 2 dependence is generated by the centrifugal part of the kinetic 
term 

(35) 

i 

Therefore if our field theory is renormalized to remove ultraviolet divergences, 
the contribution of high spin partial waves will decouple in the same manner 
as other high momentum effects such as dependence on an explicit cutoff 
mass scale. The error in neglected partial waves with spin larger than some 
value J max is roughly equivalent to the error in using a cutoff mass A eff where 

Ae ff ~%S (36) 

and t is the characteristic radius for the process being studied. 

In two-dimensional 4 theory, renormalization can be done simply by 
normal ordering the interaction. In the context of spherical field, normal 
ordering corresponds with explicitly subtracting the contribution of diagrams 
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of the type shown in Figure 3. The amplitude for this diagram is proportional 
to 

XtK lnl (fit)I H (fit). (37) 

In the Schrodinger language this corresponds with a counterterm propor- 
tional to 

XbK\ n \(iJ,t)I\ n \{nt)q-iqi. (38) 

Let us now return to the two-loop perturbative calculation. In presence 
of an explicit cutoff mass A e ff, the A e g-regulated calculation for the two-loop 
result s oc (k 2 ) includes a correction of size 0(A~g). From relation (36), we 
expect 

^(^)-^ max (^ 2 )~0(J ma 2 x ). (39) 



6 Non-perturbative self-energy 

Let us define c = and A(£) to be the difference between the two-point 
0o correlator for general c and c = 0, 



A(t) = <0| MQMty |0> c - <0| 0o(*)0o(O) |0) 



c=0 



(40) 



We can compute A(£) using (13) and (|22| ) along with normal-ordering coun- 
terterms described in the previous section. For c = 0.20 we have computed 
A(t) for J max = and J max = 1. For J max = the Schrodinger time evolution 
generator is 



i a 2 



2tdql 2 



% + ct % ~ QctK (fit)I (fit)q . 



For J n 



1 the time evolution generator is 



' 21 
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(ill) 



-Get [K (fit)I (fd) + 2K 1 {iit)I 1 {f J it)} q 2 + qft + ^ 



where 



q s i 



H+Q-i „ 
V2 ' qal 



qi-Q-i 



(41) 



(42) 



(43) 



For each value of J max we did the computation two ways, once using the 
method of finite differences and once using diffusion Monte Carlo. The results 
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are shown in Figure 4, where we have used the abbreviations fd for finite 
difference method and mc for Monte Carlo. This plot should be regarded as 
a measure of the accuracy of the Monte Carlo algorithm. The error of the 
finite difference algorithm is likely negligible as we were able to reproduce 
the two-point O correlator for free field theory to an accuracy of 0.1%. 
Let U(k 2 ) denote the (ft self-energy, defined such that the full (p propagator 

is 



f(k 2 



k 2 +fi 2 +n(k 2 ) ■ 

We have calculated -II(£; 2 ) for c = 0.05,0.10,- 
Jma.x = 2 the time evolution generator is 



i_ 

' it 
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,0.40 and J max = 2. For 
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where 



q S 2 



<?2+g-2 



Our results were calculated using diffusion Monte Carlo and are shown in 
Figure 5. The physical squared-mass /ip hys can be found by extrapolating 
— Ii(k 2 ) to its intersection point with the line k 2 + 1 (in units where the bare 
mass n = 1), 



-II( 



2 

A'phys 



ivs 



+ 1. 



One fact that is useful in performing this extrapolation is 

<0, 



dU 

' dk 2 



k 2 -- 



^phys 



(47) 



(48) 



since the residue of the pole at — ^ 2 hys must be less than or equal to l. 2 As 
Figure 5 shows (and also visible in the plots of s j max (fc 2 ) in Figure 2), the slope 
of — U(k 2 ) appears to have the wrong sign at small values of k 2 . The reason 
for this is that high spin partial waves cannot be neglected for t ~ 0(J max ), 



2 This can be proved by spectral decomposition of the two-point <f) correlator. 
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or equivalently k 2 ~ 0(J~ 2 X ). It is therefore necessary to extrapolate — II(&; 2 ) 
from values of k 2 where this breakdown does not occur. 

Extrapolating — IT(£; 2 ) through the small k 2 region is somewhat difficult 
for J max = 2, and we could achieve greater accuracy by choosing a higher 
value for J max - But, as mentioned before, our main goal here is to introduce 
to the methods of spherical field theory. A more thorough and accurate 
calculation of this and other phenomena will be done in the future. For now 
the extrapolation method we choose is simply to extend the maximum value 
of -I1(A; 2 ) to all k 2 to the left of the maximum, as shown by the horizontal 
dashed lines in Figure 5. With this simple procedure we find the following 
values of /ip hys . 



For c near 0.4 the physical mass vanishes. As discussed in the literature, [2] 
and [[|, this is associated with spontaneous breaking of </> — > —<fi reflection 
symmetry. Despite using a small value for J max and a minimal amount of 
computer labor, we find that our value for the critical coupling is in agreement 
with a recent lattice calculation, |J, which found the critical coupling at 
c = 0.407. 

7 Summary 

Spherical field theory is a new approach to studying quantum field theory. 
It treats a general (i-dimensional field theory as a system of coupled one- 
dimensional theories, which are then solved using partial differential equa- 
tions. In this paper we introduced the methods of spherical field theory using 
the example of Euclidean 4 theory in two dimensions. We first derived the 
spherical Feynman rules, showed that high spin partial waves decouple, and 
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c 




(49) 
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then, using a spin cutoff J max = 2, computed the self-energy for non- 
perturbative values of the coupling. Our results show clear evidence of a 
phase transition near c = 0.4. The methods used here are easily general- 
ized to any d-dimensional Euclidean field theory involving scalar fields. The 
extension to fermionic and gauge systems is currently being studied. 

We expect spherical field theory to have many applications in the areas of 
particle theory and statistical physics. It offers an efficient and flexible way 
to calculate non-perturbative quantities in field theory. For simple systems 
(small J max ) either the method of finite differences or Monte Carlo simulation 
can be used. For larger systems only Monte Carlo methods are feasible, but 
here again we have several qualitatively different options. In our computation 
of the 4> self-energy we used diffusion Monte Carlo to model the time evolution 
equation. Another option, though, is to compute the one-dimensional path 
integral shown in flT2p using standard lattice field methods. The ability 
to calculate quantities in several different ways is often useful in finding 
efficient algorithms and error-checking. In this regard we believe spherical 
field calculations should be useful in lattice field calculations and vice versa. 
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Figure 1. The two-loop self-energy diagram for (fi . 

Figure 2. Plot of the two- loop self-energy, with normalization given in (32), 
for J max = 0,l,2,3,oo. 

Figure 3. The only divergent diagram, which can be removed by normal 
ordering. 

Figure 4. Plots for A(t) as calculated by finite difference (fd) and Monte 

Carlo (mc) methods, with J max = 0, 1 and c = 0.20. 

Figure 5. Plots for -n(£; 2 ) with J max = 2 and c = 0.05, 0.10, • • -0.40. 
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Figure 1 
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Figure 3 
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Figure 5 
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